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SECTION  1 


INTRODUCTION 


The  accuracy  and  dispersion  of  free  flight  vehicles  has 
been  a  problem  in  aerodynamics  and  ballistics  for  many  years. 
Until  the  present  time,  the  primary  investigations  into  causes 
and  effects  of  jump  (the  angle  between  the  line  of  boresight 
and  the  line  connecting  the  point  of  launch  with  the  instan¬ 
taneous  position  on  the  trajectory)  and  dispersion  have  been 
directed  toward  projectiles  and,  in  particular,  artillery 
rounds.  A  full  program  to  investigate  jump  and  dispersion 
characteristics  of  low  trajectory  finned  bodies  has  been  lack¬ 
ing  and  therefore  is  the  subject  of  this  report.  The  purpose 
of  this  analysis  is  to  develop  a  basic  understanding  of  the 
parameters  causing  the  jump  and  dispersion  of  flechettes.  The 
flechette,  being  a  gun  launched  finned  body,  requires  a  differ¬ 
ent  approach  to  the  problem.  The  old  concept  employed  in  the 
analysis  of  the  dispersion  of  artillery  rounds  is  that  the 
dispersion  results  from  initial  launch  disturbances  imparted 
by  the  gun  to  the  shell  (References  1  and  2).  This  concept  is 
no  longer  valid  for  flechettes  since  the  flechette  is  a  fin 
missile,  sabot  launched,  and  its  dispersion  must  be  tied  to 
the  disturbances  it  encounters  when  clearing  the  muzzle  blast 
and  sabot  separation  region.  In  addition,  asymmetries  are  more 
prevalent  in  finned  bodies  than  projectiles  and  a  finned  body 
is  more  apt  to  be  influenced  by  the  blast.  These  factors  must 
be  taken  into  account  by  a  theory  involving  finned  bodies. 

In  order  to  develop  this  new  approach,  (1)  a  theoretical 
expression  for  jump  and  dispersion  had  to  be  developed,  (2) 
the  theory  had  to  be  validated,  (3)  free  flight  test  firings 
had  to  be  undertaken  and  initial  condition  data  extracted,  and 
(4)  the  test  firing  results  had  to  be  correlated  with  the  vali¬ 
dated  theory.  The  Jump  and  Dispersion  Theory  was  developed, 
in  general,  for  both  fin  and  spin  stabilized  missiles  in  air. 
The  theory  includes  the  effects  of:  initial  conditions,  magnus 
aerodynamic  asymmetries,  and  gravity.  In  the  past,  theory 
development  for  projectiles  included  only  initial  angle  of 
attack  and  initial  angular  rate  (References  1  and  3).  Initial 
transverse  velocity  was  considered  non-existent  (Reference  4)  o 
negligible.  Zaroodny  (Reference  5)  included  a  linear  momentum 
term  to  account  for  any  transverse  motion  of  the  projectile  but 
attributed  it  to  the  gun  during  recoil.  Any  transverse  impulse 
imparted  to  the  projectile  by  the  blast  was  ignored.  Other 
authors  including  Sterne  (Reference  2)  attributed  the  jump  only 
to  bore  cleararce  and  therefore  only  included,  effectively,  the 
initial  angle  of  attack.  Magnus  effects  were  always  neglected 


in  previous  studies  either  due  to  lack  of  familiarity  with  the 
subject  or  lack  of  data.  In  general,  all  cross-forces,  except 
lift,  were  neglected  mainly  for  convenience  sake.  Zaroodny, 
however,  cautioning  against  wholesale  simplifying  said  "it 
would  seem  desirable  that  our  formulas  allow  us  to  include 
these  other  forces  as  the  experimental  information  on  these 
forces  becomes  available."  Aerodynamic  asymmetries  were 
neglected  for  projectiles  but  included  in  Murphy’s  work  (Refer¬ 
ence  6).  It  was  not  until  Nicolaides  (References  7,  8, and  9) 
that  all  four  factors  affecting  dispersion;  initial  angle  of 
attack,  initial  angular  rate,  initial  transverse  position  and 
velocity,  were  put  into  one  theory.  The  work  presentee  here 
expands  the  work  of  Nicolaides  to  include  all  parameters  affect¬ 
ing  dispersion  in  detail.  Three  separate  equations  comprise 
the  theory  to  include  the  complete  range  of  roll  rates.  Before, 
only  high  roll  rates  were  considered;  with  the  study  of  finned 
bodies,  the  roll  rate  range  extends  down  to  zero  roll  and  accu¬ 
rate  theories  had  to  be  deduced  from  known  aerodynamic  equations. 

To  validate  the  theory,  a  s’ x-degree-of- freedom  trajectory 
computer  program  numerically  integrating  the  equations  of 
motion  was  utilized  (References  10,  11,  and  12).  The  valida¬ 
tion  consisted  of  four  phases.  The  procedure  began  with  the 
most  basic  theory  equation  and  consecutively  added  terms  tc 
validate  the  entire  theory.  Initial  conditions,  magnus ,  asym¬ 
metries  and  gravity  were  successively  validated  with  roll  rate 
and  velocity  varied  in  each  phase. 

Before  the  advent  of  adequate  photographic  material,  obtain¬ 
ing  test  data  was  often  difficult.  At  first,  jump  target  data 
was  taken  separately  from  yaw  data.  The  thinking  was  that  the 
yaw  data  was  part  of  the  projectile's  characteristics  and  not 
affecting  jump.  As  photographic  methods  improved,  and  theories 
developed,  the  data  was  correlated.  The  correlation  of  the 
data  was  often  a  problem.  A  fit  of  the  motion  to  a  least 
squares  method  was  difficult.  Fowler,  Kent,  and  Hitchcock 
developed  a  method  that  would  plot  the  magnitude  of  the  yaw 
separately  from  the  orientation  and  then  fit  the  curves  sepa¬ 
rately.  A  better  method  was  developed  by  McShane-Charters  Turet¬ 
sky  that  approximated  the  yawing  motion  to  a  circle.  For  pro¬ 
jectiles  the  method  has  been  refined  and  is  an  excellent  method. 
However,  for  finned  bodies  with  not  always  circular  angular 
motions,  a  different  method  of  data  analysis  had  to  be  devised. 
Utilizing  the  free  flight  data  taken  by  test  engineers  at 
Frankford  Arsenal  on  a  number  of  flechettes,  the  least  squares 
method  was  employed  to  fit  the  data  presented  here.  The  nearly 
planar  oscillations  of  the  flechette  in  the  first  few  feet 
downrange  were  fit  to  a  pure  pitching  motion  (References  13 
and  14)  and  the  position  downrange  fit  to  a  third  order 
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polynomial.  From  these  results,  angle  of  attack,  angular  rate 
and  transverse  position  and  velocity  were  determined  for  the 
first  few  feet  downrange.  Before,  there  was  some  controversy 
as  to  whether  or  not  the  least  squares  fit  could  be  extrapo¬ 
lated  back  to  the  muzzle.  Zaroodny  contended  that  the  x=0 
position  had  to  be  taken  out  of  the  blast  region  to  allow  the 
aerodynamic  equations  to  be  valid.  On  the  other  hand,  Kent, 
Hitchcock,  Fowler  and  Sterne  held  to  the  fact  that  the  free 
flight  region  began  the  instant  the  projectile  left  the  bore. 

In  the  analysis  of  flechettes  the  position  x*-0  is  taken  some¬ 
where  downrange  after  the  sabot  separation  sequence  has  occurred. 
This  is  seen  to  be  3  to  5  feet  downrange  and  assumed  clear  of 
any  muzzle  blast  effects. 

The  striking  shortcoming  of  previous  works  is  the  lack  of 
correlation  between  test  data  and  valid  theory.  For  the  flech- 
ette,  correlation  between  the  theory  and  test  data  was  under¬ 
taken  as  well  as  correlation  between  test  results  and  first 
maximum  yaw  data.  Currently,  the  first  maximum  yaw  theory 
(Reference  15)  is  held  by  some  to  be  an  accurate  method  of 
predicting  dispersion.  This  theory  disallows  any  influence  of 
initial  angular  rate,  transverse  position,  or  velocity  on  dis¬ 
persion.  The  dispersion  analysis  present  1  here  disproves 
this  theory  with  actual  test  data.  The  deails  of  each  of 
these  aspects  of  this  program  are  developed  in  the  following 
sections . 
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SECTION  II 


DISPERSION  THEORY 


Dispersion  relationships  for  free  flight  vehicles  are 
embedded  in  the  trajectory  equation  of  any  such  aeroballistic 
body.  To  evaluate  the  trajectory  equation  and  thus  the  disper¬ 
sion,  the  linear  second-order  differential  equation  of  angular 
motion  is  a  logical  starting  point. 


w  +  N,w  +  N0 w  =  N  e ipt  +  Ff, 
12  3  4 

where  N^,  N^,  and  are  constants. 


(1) 


In  this  discussion  of  dispersion  theory,  it  is  assumed  that, 

1. 

theory  development. 


total  velocity,  u  ,  is  constant,  equal  to  u  in  the 


2.  all  force  and  moment  coefficients  dependent  on 
angle  of  attack  are  considered  to  be  linear  with 
angle  of  attack. 


3. 


all  force  and  moment  coefficients  independent  of 
angle  of  attack  are  considered  to  be  constant. 
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4.  a  linear  relationship  exists  between  x  (distance 
downrange)  and  time  for  the  nondrag  case. 


5.  roll  rate,  p,  is  considered  to  be  constant. 

6.  products  of  force  and  moment  derivatives  are  negli' 
gible,  except  those  involving  and  . 

e  e 

Utilizing  these  assumptions,  and  the  jinomial  expansion  of 
(Z^-m)’1,  2,  3,  4,  and  5  become: 
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The  solution  to  Equation  1  is  that  of  tricyclic  motion;  that  is, 
w  =  K  e^l1  +  K2e<l>2t  +  K3e  iPt  +  K4 
where  the  compler.  coefficients  are; 

-  _  w0  ’(^2,1)  wg  +  K3  ((t>2.riP) 

Ni  2  "  - - 

^1,2  ‘^2,1 


K  =  N3 
3  "  (ip-<>  1)  (ip-02) 
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k  -J5l 
4  N2 


(9) 


and 


(10) 


The  trajectory  equation  for  free -flight  motion: 

S  =  (\£  *  iuq) 

An  expression  for  q  is  obtained  from  the  equations  of  motion 


(11) 
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yielding  a  solution  of  the  form: 


ipt  o 
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+  k3e  +  k4t  +  k5t  +  k6 


(13) 


where  the  entire  expression  for  the  solution  is: 
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HIGH  ROLL  RATE  THEORY 


For  roll  rates  greater  than  100  rad/sec,  Equation  16  reduces  to 
an  approximate  solution.  Integration  of  the  double  integral  gives: 

.**  .  _t_ .  i 

o- o'  (ip/  ip  (ip? 

For  high  roll  rates,  the  first  and  third  terms  go  to  zero,  leaving  only 
the  second  term  to  affect  dispersion.  Applying  this  approximation  to 

i 

Equation  16  : 


t  t 

J  =  ” 


(17) 
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(18) 


where,  by  applying  previous  aerodynamic  relationships: 


(19) 
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Substituting  19  and  20  into  18  and  expanding  the  various  terms: 


The  mil- relat ion  offers  a  method  to  define  the  Jump  Angle  from 
Equation  22. 


Jump  Angle  *  ~  (10^) 
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Equation  24  gives  an  approximation  for  the  Jump  Angle  for  high 
roll  rate  cases  with  gravity,  at  any  position  x  down  range. 

LOW  ROLL  RATE  THEORY 

For  roll  rates  less  than  100  rad/sec  but  having  a  parameter,  pt, 
greater  than  1 ,  Equation.  16  can  be  reduced  to  another  approximation. 
As  before,  integration  of  the  double  integral  yields  Equation  17 


n 


OP ? 


For  low  roll  rates  all  three  terms  are  significant  to  dispersion. 
Equation  16  now  becomes: 

/^w+*P^na  T*-  ZifiA/  1  fK  4 /Zw+ ipZno\  lor"!  9 


r3  ♦  MVi-eiPt)  j,  +  f-}2 

^  ^  ^  (25) 

£i ,  «A  ii ,  Mil  t 

m  A*i  *2/  ip  mp  J 

The  arm,  or  rolling  trim  vector  must  be  separately  examined. 
From  Equation  8,  10 
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Numerical  inspection  of  the  three  denominator  terms  indi¬ 
cates  that  the  first  two  terms  can  be  neglected.  Each  term  is 
not  only  less  than  1  percent  of  the  third  term  but  also  they're 
subtracted  from  one  another  to  make  their  contribution  even  more 
minimal.  Thus  is  approximated  by, 

i  [ ipCv*  0-  Cm^  1 

md 


(S+  -mTd(£)c*Pe]+i  -&Y 


for  low  roll  rates,  the  second  term  in  the  numerator  and  the 
first  term  in  the  denominator  dominate  all  other  teinis  and 
become  the  only  significant  terms.  Thus, 


-  ulCMi,‘. 

K3“-  -CtJ— 


~i *ie  same  approximation  holds  true  for  applicable  terms  in 
Equation  25  ,  thus  reducing  the  jump  angle  equation  to: 
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Combining  terms  and  dropping  the  negligible  second  last  term, 


Equation  28  accurately  approximates  the  jump  angle  for  roll  rates: 


p  <  100  rad/sec 

pt>  1.0 

VERY  SLOW  ROLL  RATE  THEORY 

For  very  low  roll  rates;  that  is,  p  >  0  and  pt  <  l,  Equa¬ 
tion  28  is  again  applicable.  ~ 

One  approximation  is  used,  however,  and  that  is  that  cos 
and  sin  are  approximated  by  power  series. 


Substituting  and  simplifying 


SECTION  III 


VALIDATION  OF  THEORY 


The  theoretical  expressions  for  Jump  Angle;  Equations  24, 
28,  and  30;  show  that  the  dispersion  depends  on  the  initial  con¬ 
ditions,  aerodynamic  coefficients,  distance  downrango ,  and  mass 
parameters.  Dispersion  for  this  theoretical  analysis  is  defined 
to  be  the  deviation  from  the  line  of  fire.  By  analyzing  only 
one  flechette  configuration  to  validate  the  theory,  the  produc- 
ibility  Ground  Point,  and  taking  all  cases  to  be  evaluated  at 
1000  feet  downrange ,  then  the  expression  for  the  Jump  Angle  can 
only  be  affected  by  the  initial  conditions  and  aerodynamic  coef¬ 
ficients  . 

To  assure  that  the  three  equations  for  Jump  Angle  are  valid 
and  to  show  the  effects  for  various  initial  conditions  and  aero¬ 
dynamic  coefficients,  the  expressions  for  the  Jump  Angle  were 
evaluated  for  a  series  of  cases  and  compared  to  numerical  inte¬ 
gration  of  the  six-degree-of- freedom  equations  of  motion,  (6-D) . 
The  series  of  cases  is  broken  down  into  various  phases  of  devel¬ 
opment.  Phase  I  considers  various  initial  conditions  but  with 
only  the  restoring  and  damping  aerodynamic  coefficients.  This 
phase  validates  the  use  of  initial  conditions  alone.  Phase  II 
utilizes  a  set  of  constant  initial  conditions,  except  for  roll 
rate,  and  constant  restoring  and  damping  coefficients,  while 
varying  magnus  coefficients  to  determine  their  influence.  Phase 
III  brings  into  consideration  all  che  aerodynamic  coefficients 
to  include  the  configurational  asymmetry  coefficients.  Differ¬ 
ent  coefficients  are  used  by  varying  the  initial  velocity  and 
roll  rates  are  varied  to  evaluate  high,  low,  and  very  low  roll 
theories.  Phase  IV  considers  the  effects  of  gravity  for  various 
initial  velocities  and  roll  rates.  No  configurational  asymmet¬ 
ries  are  used  in  order  to  isolate  the  gravitational  influence. 
Values  for  all  coefficients  are  found  in  the  Appendix  as  well 
as  other  data  including  mass  parameters.  Since  computations 
were  done  at  1000  feet  downrange,  the  Jump  Angle  in  mils  is 
equivalent  to  the  deviation  from  the  line  of  fire  in  feet  for 
all  presented  cases.  The  axis  system  used  throughout  this  anal¬ 
ysis  is  illustrated  in  the  list  of  symbols. 

PHASE  I 

To  validate  the  effects  of  initial  conditions  with  restor¬ 
ing  and  damping  coefficients  only,  36  cases  were  evaluated  using 
the  high  roll  rate  theory,  Equation  24.  The  cases  are  divided 
into  four  sections  isolating  different  initial  conditions  and 
their  effects. 
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Cases  1-9 


The  first  section  shows  the  effects  of  roll  rate  and  veloc- 

-t*-  _ t  -r* 

ity  with  zero  S  .  a  and  a. 

Table  1  clearly  indicates  that  no  deviation  from  the  line 

of  fire  occurs  if  "s’  ,  ""a*  ,  and  "oT  are  set  to  zero.  Roll  rate 

oo’  o 

and  velocity  changes  have  no  effect  on  the  Jump  Angle  for  this 
particular  situation.  This  is  a  trivial  solution,  it  being 
obvious  from  inspection  of  Equation  24. 

TABLE  1.  THEORY  VALIDATION,  RESTORING  AND 
AND  DAMPING  MOMENTS,  CASES  1-9 


c 

A 

Coefficients 

Initial  Conditions 

CyE 

CZE 

CME 

CNE 

J.  A  (mils) 

S 

E 

K 

a 

o 

K 

Po 

uo 

CM  a 

Cm  +  cM 

4  a 

CZpg 

CMpp 

6-D 

Theory 

1 

9 

0 

9 

31416 

■ 

i 

0  +  01 

0  + 

Oi 

2 

0 

0 

0 

18850 

5000 

0  +  Oi 

0  + 

Oi 

j 

0 

0 

0 

6283 

0  +  Oi 

0  + 

Oi 

4 

0 

0 

0 

31416 

0  +  Oi 

0  + 

Oi 

5 

•0 

0 

0 

18850 

3000 

A1 

c 

D 

0  +  Oi 

0  + 

Oi 

6 

0 

0 

0 

6283 

0  +  Oi 

0  + 

Oi 

.7 

0 

0 

0 

31416 

0  +  Oi 

0  4. 

Oi 

8 

0 

0 

0 

18850 

1000 

0  +  Oi 

0  + 

Oi 

9 

0 

0 

0 

6283 

* 

_ ! 

0  +  Oi 

0  + 

Oi 

Cases  10-18 

The  second  section  gives  the  effects  of  initial  trans¬ 
lational  velocity,  SQ  =  y  +  iz,  with  various  roll  rates  and 

velocities.  To  assure  the  solution  is  correct  in  three  dimen¬ 
sional  space,  the  initial  translation  velocity  is  given  in  both 
y  and  iz  directions.  Equation  24  reduces  to: 
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J .  A. 


1000  T* 
u  bo 


The  correlation  between  the  theory  and  the  6-D  integration 
for  Cases  10-18  is  excellent  as  shown  in  Table  2.  The  Jump 
Angle  is  seen  to  be  affected  by  velocity  but  not  roll  rate,  as 
would  be  expected  from  the  reduced  Jump  Angle  equation.  Figure 
1  illustrates  the  deviation  from  the  line  of  five  for  initial 


TABLE  2.  THEORY  VALIDATION,  RESTORING  AND 
DAMPING  MOMENTS,  CASES  10-18 


J. A.  (mils) 

6-D 

Theory 

20.002  + 
20.0061 

20.000  + 
20.000  1 

20.002  + 
20.0061 

20.000  + 
20.000  i 

20.002  + 
20.0061 

20.000  + 
20.000  1 

33.346  + 
33. 368  i 

33.333  + 
33. 333  i 

33.346  + 
33.3681 

33.333  + 
33.333  i 

33.346  + 
33. 368  i 

33.333  + 
33.333  i 

100.254  + 
100. 765  i 

100.000  + 
100.000  i 

100.254+ 
100. 764 i 

100.000  + 
100.000  1 

100.257  + 
100.7651 

100.000  + 
100.000  i 

c 

A 

S 

E 


10 


11 


12 


13 

14 


15 


16 


17 


18 


Initial  Conditions 


Qtr 


100+ 

lOOi 


100  M 
100  i 


100+] 
100  i 


100+1 
100  i 


100+1 
100  i 


100+1 
100  i 


100+] 
100  i 


100+1 
100  i 


100+1 
100  i 


0 


0 


0 


0 


0 


0 


0 


0 


Otr 


0 


0 


0 


0 


0 


0 


0 


0 


o 


31416 


18850 


6283 


31416 


18850 


6283 


31416 


18850 


6283 


LL 


5000 


3000 


1000 


Coefficients 


CZa 
cMa 

cMq+CMA 


A1 


'P0 


Cm 


P3 


o 


CzE 

CMH 

CNE| 


0 


16 


Cases  19-27 


The  third  section  gives  the  effects  of  initial  angle  of 

attack,  aQ,  with  various  roll  rates  and  velocities.  Again  a 

complex  initial  condition  is  used  to  validate  the  theory  in 
three  dimensional  space.  Equation  24  reduces  to: 


J .  A.  = 


°\ai/ 


CM  -  i(plx)  C7 

Ma  Vi USj  2  a 


1000 


Table  3  shows  the  range  of  error  between  the  6-D  computa¬ 
tion  and  the  theory  to  be  0.036  to  0.040  mils  in  the  y-direction 
and  0.038  to  0.041  mils  in  the  z-direction.  Although  the  y- 
direction  deviations  differ  in  sign,  the  error  between  them  is 
approximately  0.00225  degrees,  an  extremely  small  angle.  This 
angle  will  give  an  approximate  deviation  of  0.04  feet  from  the 
line  of  fire  at  1000  feet  downrange.  With  the  J.A.  being  so 
close  to  zero  it  can  be  expected  that  the  signs  may  differ  due 
to  computational  errors.  The  results  do  show  Jump  Angle  vari¬ 
ance  with  both  roll  rate  and  velocity.  The  largest  changes 
occur  as  velocity  goes  to  1000  ft/sec. 

Cases  28-36 

The  fourth  section  gives  the  effects  of  initial  angular 

m  fci 

rate  dQ,  with  varying  roll  rate  and  velocity.  An  angular  rate 

of  250  rad/sec  is  used  in  both  directions  of  the  complex  plane 
to  test  validity  in  three  dimensional  space.  Equation  reduces 
to : 


J.A.  = 


“»  (ka) 


1000 


Table  4  indicates  excellent  agreement  between  the  theory 
and  6-D  computations.  Roll  rate  is  found  not.  to  affect  the 
Jump  Angle  appreciably  but  velocity  does,  as  would  be  expected 
from  the  reduced  Jump  Angle  equation.  Figure  3  shows  the  dis¬ 
persion  pattern  while  Figure  4  illustrates  the  trajectories. 
Cases  28,  29,  and  30  are  plotted  as  one  point  due  to  the  small 
difference  between  the.  Cases  31,  32,  33,  and  34,  35,  and  36 
are  plotted  similarly. 


TABLE  4 .  THEORY 


DAMPIN 


c 

A 

S 

E 

Initial  Conditions 

9 

S0 

a 

Po 

uo 

28 

0 

0 

250+ 

250i 

31416 

.5000 

29 

0 

0 

2^0+ 

250i 

18850 

30 

0 

0 

250+ 

250i 

6283 

31 

0 

0 

250+ 

250i 

314  3  6 

3000 

32 

0 

0 

250+ 

250i, 

18850 

33 

0 

0 

250+ 

250! 

6283 

34 

0 

0 

250+ 

250i 

31416 

1000 

35 

0 

0 

250+ 

250i 

18850 

36 

0 

. 

0 

250+ 

250i 

6283 

VALIDATION,  RESTORING  AND 
i  MOMENTS,  CASES  28-36 


J.A.  (mils) 

6-D 

Theory 

i 

-2.027 

-2.034i 

-2.073 
-2. 073! 

-2.025 

-2.030i 

-2.073 
-2. 073i 

-2.027 
-2. 029i 

-2.073 
-2. 073i 

-1.961 

-1.9671 

-1.970 
- 1 . 970i 

-1.962 

-1.9661 

-1.970 
- 1 . 970i 

-1.964 

-1.9641 

-1.970 

-1.970! 

-5.238 
-5, 274i 

-5.540 
-5. 540i 

-5.243 
-5  264 i 

•  5 . o40 
-5.5401 

-5.254 

-5.2601 

-5.540 

-5.540i 

Coefficients 


CZ 

Cm 


a 

a 


CMq+°M^ 


CZ 


P0 


Cm 


P0 


We 

CZE 

CME 

CNE 


A1 


0 


0 
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CASES  28,29. 30 
CASES  31,32,33 
CASES  34,35,36 


PHASE  II 


To  validate  the  effect  of  Magnus  Forces  and  Moments  on  the 
dispersion  of  flechettes,  21  Cases  were  run  varying  the  initial 
roll  rate  and  magnus  coefficients.  All  other  conditions  were 
held  constant.  The  variance  of  magnus  coefficients  with  Mach 
number  had  to  be  chosen  since  no  data  was  available.  Arbitrar¬ 
ily,  the  ratio  of  C_  /CM  was  chosen  to  be  the  same  as  that 

2P0  MP0 

of  C  /Cw  .  The  magnus  coefficients  used  are  presented  as 
a  a 

functions  of  Mach  Number  in  the  Appendix  with  only  the  values 
at  Mach  4.5  tabulated  here  for  identification  sake: 

TABLE  5.  MAGNUS  COEFFICIENTS,  AT  MACH  4.5 


t7 

P0 

+  34.8 

+  110.0 

+  31.6 

*  100.0 

+  28.4 

+  90.0 

Equation  24  now  becomes: 


Initial  conditions  used  in  this  section  are  consistent  with 
those  of  other  sections  to  provide  a  basis  for  comparison. 
Three  cases  of  zero  magnus  were  run,  one  at  each  roll  rate  to 
provide  a  standard  to  judge  the  influence  of  magnus. 


24 


iff  . 


The  effects  of  magnus  coefficients  on  dispersion  are  mini¬ 
mal  as  seen  in  Table  6.  The  variance  between  the  zero  magnus 
cases  and  any  other  case  is  found  not  to  be  greater  than  0.209 
mils  (or  feet  at  1000  feet  of  range) .  *n  order  to  obtain  the 
maximum  magnus  effects,  the  largest  possible  magnus  coefficients 
were  used.  Hence,  C7  =  34.8  and  CM  ~  J  10.0  are  the  largest 

^pB  Mpb 

possible  coefficients  since  cases  40  and  *9  become  unstable. 

Table  6  indicates  the  effects  (for  positive  magnus  coefficients). 

1.  increasing  horizontal  dispersion  with  increasing  p 

2.  decreasing  vertical  dispersion  with  increasing  p 

3.  increasing  horizontal  dispersion  with  increasing 
magnus 

4.  decreasing  vertical  dispersion  with  increasing 
magnus 

(for  negative  magnus  coefficients) 

5.  decreasing  horizontal  dispersion  with  increasing  p 

6.  increasing  vertical  dispersion  with  increasing  p 

7.  decreasing  horizontal  dispersion  with  decreasing 
magnus 

8.  increasing  vertical  dispersion  with  decreasing 
magnus 

For  example,  Figure  5  illustrates  the  effects  of  roll  rate  for 
constant  magnus  coefficients  of  +_  90^  (1,2, 5, 6  above).  Figure 
6  illustrates  the  effects  of  magnus  for  a  constant  sample  roll 
rate  (3, 4, 7, 8  above).  Obviously,  when  only  a  0.209  mil  maxi¬ 
mum  deviation  due  to  magnus  occurs  when  the  situation  is  geared 
toward  finding  the  largest  effect  due  to  magnus,  smaller  devia¬ 
tions  due  to  magnus  would  be  found  in  actual  situations.  It 
can  be  concluded  that  Magnus  has  no  large  effect  on  dispersion 
although  it  could  be  significant  if  the  total  dispersion  is 
close  to  zero. 

PHASE  III 

To  validate  the  effects  of  aerodynamic  asymmetries  on  dis¬ 
persion  of  flechettes,  a  large  number  of  cases  were  run  varying 
roll  rate,  velocity,  and  initial  conditions  while  holding  the 
asymmetry  coefficients  constant.  The  asymmetries  coefficients 
were  selected  to  allow  1  degree  of  non-rolling  trim  to  exist 
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37  Cz  =  0.0 

,38  cM  =0.0 

1 39 


40  Cz  =  34.8 

41  „  ^ 

cMn  O  =  nao 


C~  =  3L6 
ZP{3 

Cm  «  100.0 
45  pe 


CZp3  =  28.4 
Cm  *  90.0 


49  Cz  =-34.8 

50  ^ 

51  CMp0-nO.O 


52  CZ  =-  31.6 

PP 

53  cm  =-100.0 

54  P3 


55  Cz  =-  28.4 

56  Pe 

Cm  =  -  90.0 

57  P0 


31416  rad/sec  18850  rad/sec 


(6-D)  y 

17.994+ 

18.042L»Theory 
17.900+ 
17.954 i 


Unstable 


18.203+ 

17.8491 


17.899 

17.9541 


18.183+ 

17.8691 


17.899 

17.954i 


Unstable 


17.807+ 

18.2581 


17.899+ 
1 7 . 954 i 


17.826+ 
18. 2331 
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17.9541 
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18. 1701 
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17.9421 
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17.090+ 
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17.9791 
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17.9311 
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17.9871 
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17.977+ 
18. 059i 


17.920+ 

17.9311 


17.920+ 

17.9311 


17.920+ 

17.9311 
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Figure  5,  Dispersion;  Phase  IT  Cases 
46,  47,  48.  58.  56,  and  57 


while  the  flechette  was  in  flight.  The  asymmetry  coefficients, 
CyE>  C^g,  Cj^g  and  C^g  are  presented  in  the  Appendix  as  a  func¬ 
tion  of  mach  number.  The  variance  with  mach  number  was  chosen 
arbitrarily:  the  ratio  of  asymmetry  force  to  asymmetry  moment 
identical  to  the  ratio  of  to  .  The  wide  range  of  roll 

rates  makes  mandatory  use  of  all  three  dispersion  theories. 

The  governing  equations  are  presented  as  they  apply. 


Cases  58-90 


The  first  set  of  cases  utilizes  zero  initial  disturbances 
while  varying  velocity  and  roll  rate.  For  roll  rates  of  31416 
rad/sec  down  to  100  rad/sec  the  High  Roll  Rate  Theory  yields 
the  governing  equation, 


8  m 


1000 


For  roll  rates:  p<  100  rad/sec  and  pt  >  1.0,  the  Low  Roll  Rate  Theory 
takes  effect: 


Finally,  the  very  Slow  Roll  Rate  Theory  applies  for  values  of  pt<  1.0: 


Tables  7,  8,  and  9  list  Cases  58-90: 
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TABLE  7.  THEORY  VALIDATION,  ASYMMETRIES, 
CASES  S  8  -  6  8 


Coefficients 


Initial  Conditions 


a  a 
o  o 


°za 

CM„ 


CM  +  CM 


P*3  C* 


CYE 

CZE 


P£  CNE 


0  0  0  314'6 


0  00  18850 


0  00  6283 


0  0  0  500 


000  300 

0  0  0  100  5000 


0  0  0 


0  0  0 


0  0  0 


0)0  0 
0  !  0  0 


J. A.  (mils) 


6-D  Theory 


0.018- 

0.0131 

0.030- 
0.027i 
0.060-< 
0.  1271 

0.997- 
0. 992i 
1.620- 
1.7211 
4.574- 
4. 896i 
8.666- 
12. 280i 
20.669- 
26.4 18i 

-7.973 
-62. 1971 


0.018- 

0.0  Ui 

0.029- 
0. 025i 

0.064- 
0.  1301 
1.013- 
1.0091 
1.688- 
1.6831 
4. 675- 
4. 9751 
8.780- 
12.4891, 
’  21.  150- 
26. 927i 

-8.210 
-63. 2  lOi 


-29.857  -30.372 
-61.459i  -62,3531 
-49.706  -50.427 
-49. 706i  -50. 4271 


ft* 


TABLE  8.  THEORY  VALIDATION,  ASYMMETRIES 
CASES  69-79 


Initial  Conditions 


0  31416 


0  18850 


0  0 


6283 


Coefficients 


CZa  CYE 

■v0"*  4 


100  3000 
50 


J.A.  (mils) 


6-D 

Theory 

0.008- 
0. 004i 

0.009- 

0.0041 

0.013- 

0.008i 

0.013- 
0. 008i 

0.033- 
0. 028i  1 

0.034- 

0.029i 

0.394- 
0. 3981 

0.401- 
0.  3961 

0.663- 
0.  659i 

0.666- 

0.6621 

1.841- 

1.9931 

1.994- 

1.9891 

3.780- 
4. 5131 

3.411- 
4. 1641 

5.676- 
8. 457i 

5.721- 
8. 5 16i 

9.203- 

32.6281 

9.217 

32.8971 

-10.029- 

-41.9851 

-42.2731 

-42.2731 

-33.014  -33.  1*34 
-33.014i-33.  1941 


TABLE  9.  THEORY  VALIDATION.  ASYMMETRIES 
CASES  80-90 


Initial  Conditions 


§  <T  OL. 

o  o  o 


80  0  0  0  31416 


81  0  0  0  18850 


82  0  0  0  6283 


83  0  C  0  500 


84  0  0  0  300 


85  0  0  0 


86  0  0  0  50 


87  0  0  0  25 


88  0  0  0  10 


89  0  0  0 


90  0  0  0 


Coefficients 


1  CZ0 

CM  a 

CzP3 

CMq+cM^ 

Cmpp 

J. A.  (mils) 


Theory 


Unstable 


Unstable 


0.023- 

0.0141 


1.177-  1.174- 

1. 160i  1.1651 


4.699- 
4. 702i 


10.224-  10.177- 
14. 447i  14. 476i 


24.402-  24.516- 
31 . 0 1 3i  31.2121 


-58.711  -58.450 
-58.7111  -58.450i 


Evident  from  Tables  7,  8,  and  9  is  the  fact  that  roll 
rate  has  tremendous  influence  on  the  dispersion  of  flechettes 
with  aerodynamic  asymmetries.  Figures  7,  8.  and  9  illustrate 
the  dispersion  pattern  for  these  cases.  The  6-D  computations 
and  theory  are  in  very  good  agreement  considering  the  large 
deviations  involved.  It  should  be  noted  that  the  actual  fle- 
chette  with  its  velocity  approaching  5000  ft/sec  is  affected 
very  little  by  aerodynamic  asymmetries.  However,  if  the  fle- 
chette  were  only  to  roll  very  slowly,  large  dispersion  ranges 
in  excess  of  60  mils  could  occur.  Velocity  also  has  a  notice¬ 
able  effect  on  dispersion.  Figure  10  shows  the  three  theory 
curves  from  Figures  7,  8,  9  in  composite  to  illustrate  veloc¬ 
ity  effects.  A  sample  trajectory,  Case  79,  is  shown  in  Figure 
11,  illustrating  the  curved  path  of  flight.  This  is  typical 
of  trajectories  involving  aerodynamic  asymmetries. 

Cases  91-123 

To  show  the  relation  between  the  effects  on  dispersion 
for  initial  transverse  velocity  and  aerodynamic  asymmetries  a 
second  set  of  cases  were  run.  Roll  rate  and  velocity  were 

varied  as  in  the  first  set  of  cases, ^but  SQ  was  set  at 

(100  =  lOOi)  ft/sec  with  ~aQ  -  0  and  cto  =  0.  Tables  10,  11, 

and  12  list  the  results.  For  high  roll  rate  cases,  Equation 
24  becomes: 


For  low  rate  cases,  Equation  28  becomes: 
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Comparing  Cases  91,  92,  93  in  Table  10  with  Cases  10,  11,  12  in 
Table  2  and  Cases  58,  59,  60  in  Table  7  it  can  be  concluded  that;  except 
for  possible  computational  error,  Cases  91,  92,  and  93  are  the  algebraic 
sum  of  Cases  10,  11,  12  and  58,  59,  60;  that  is,  for  example,  Case  91 
equals  Case  10  plus  Case  58.  This  fact  is  obviously  true  of  the  theory 
equations  and  is  here  shown  to  be  the  case  for  the  6-D  computations  as 
well.  Similar  comparisons  can  be  made  with  corresponding  cases  in  Tables 
2,  8,  11,  and  2,  9,  12.  Thus,  the  effects  of  aerodynamic  asymmetries  and 
those  of  initial  transverse  velocity  are  independent  of  one  another. 

Figures  12,  13  and  14  illustrate  the  Cases  91-123.  The  curves  are  of 
.the  same  form  as  Figures  7,  8,  and  9  but  differ  with  the  addition  of 

S  i  .  Maximum  effect  of  all  parameters  is  desired.  Cases  113,  114  and  115 
show  the  limit  of  parameter  combinations  by  113  and  114  going  unstable. 

Roll  rate  effects  are  again  large  and  velocity  effects  are  larger  than  in 
Cases  58-90.  Figure  15  shows  this  to  be  true  and  also  shows  the  cases 
involving  U  *  3000  ft/sec  to  be  ones  of  smallest  dispersion.  Such  was  the 
case  in  Figure  10.  Figure  16  illustrates  a  sample  trajectory,  Case  101. 

Cases  124-156 

To  establish  the  relationship  between  the  effects  on  disj>ersion  for 
aerodynamic  asymmetries  and  initial  angle  of  attack,  a  third  set  of  cases 
were  run.  Again  roll  rate  and  velocity  were  varied  as  done  previously  but 

aQ  was  set  at  (1+i)  degrees  with  SQ-0  and  ao=0.  Tables  13,  14,  and  15  tabu¬ 
late  the  results.  For  all  high  roll  rate  cases,  Equation  24  reduces  to: 
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For  very  low  roll  rates,  Equation  30  reduces  to: 


Only  for  high  roll  rates  does  the  aQ  term  appear.  aQ  should  have  no 
noticeable  effect  on  dispersion  for  p  <  100  rad/sec. 

Comparing  Cases  124,  125,  126  in  Table  13  with  Cases  19,  20,  21  in 
Table  3  and  Cases  58,  59,  60  in  Table  7  it  can  be  concluded  that  Cases 
124,  125  and  126  are  the  algebraic  sum  of  Cases  19,  20,  21  and  58,  59.  60; 
that  is,  for  example,  Case  124  equals  Case  19  plus  Case  58.  This  is 
obvious  from  the  reduced  theoretical  equations  for  Cases  124-156.  It  is 
shown  here  to  be  also  true  for  the  6-D  computations;  allowing  for  some 
computational  error.  Similar  comparisons  can  be  made  with  corresponding 
cases  in  Tables  3,  8,  14  and  3,  9,  15.  Thus  the  effects  of  aerodynamic 
asymmetries  and  those  of  initial  angle  of  attack  are  independent  of  one 
another. 

Figures  17,  18  and  19  illustrate  Cases  124-156.  The  curves  are  very 
similar  to  those  in  Figures  7,  8  and  9  with  the  only  difference  being  the 
very  small  a  contribution  in  Figures  17,  18  and  19.  Cases  135,  146,  and 
147  result  in  instabilities,  indicating  that  maximum  effect  of  the  various 
parameters  has  been  accomplished.  Effects  of  roll  rate  are  essentially 
the  same  as  in  Case  58-90  and  effects  of  velocity,  Figure  20,  the  same  as 
in  Figure  10.  Cases  with  U  -  3000  ft/sec  again  have  the  smallest  disper¬ 
sion.  Figure  21  shows  a  typical  trajectory,  Case  154. 

Cases  157-189 


To  validate  the  relationship  between  the  effects  on  dispersion  for 
aerodynamic  asymmetries  and  those  of  initial  angular  rate,  a  fourth  set 
of  cases  were  run.  As  before,  roll  rate  and  velocity  were  varied,  but 

to  set  at  (250  +  250i)  rad/sec  with  SQ  =  0  and  “a  =  0.  Tables  16,  17,  18 

give  the  results.  For  high  roLl  rates,  the  governing  equation  becomes: 
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For  low  roll  rates,  the  governing  equation: 


For  very  slow  roll  rate,  the  governing  equation: 


Comparing  Cases  157,  158,  and  159  in  Table  16  with  Cases  28,  29,  30  in 
Table  4  and  Cases  58,  59,  60  in  Table  7,  it  can  be  concluded  that  Cases 
157,  158,  and  159  are  the  algebraic  sum  of  Cases  28,  29,  30  and  58,  59, 

60;  that  is,  for  example,  Case  157  equals  Case  28  plus  Case  58.  This  is 
obvious  from  the  reduced  theoretical  equations  for  Cases  157-189.  Here 
it  is  shown  to  be  true  for  6-D  computations  also.  Any  discrepancy  can  be 
attributed  to  computational  error.  Similar  comparisons  can  be  made  with 
corresponding  cases  in  Table  4,  8,  and  17.  Thus  the  effects  of  aerodyna¬ 
mic  asymmetries  and  those  cf  initial  angular  rate  are  independent  of  one 
another. 
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33,  3251 
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34.000+ 
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31.3611 
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31.3441 
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28. 850i 

36.744+ 
29. 1691 

39.029+ 

24.8921 
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42.575+ 
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Figure  13.  Dispersion;  Phase  III 
Cases  102-112 
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TABLE  13.  THEORY  VALIDATION,  ASYMMETRIES , 
CASES  124-134 


TABLE  14.  THEORY  VALIDATION,  ASYMMETRIES, 
CASES  135-145 
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FIGURE  17.  DISPERSION:  PHASE  III 
CASES  124-134 
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Figure  18.  Dispersion:  Phase  III 
Cases  135-145 
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TABLE  17.  THEORY  VALIDATION,  ASYMMETRIES, 
CASES  168-178 


Initial  Conditions 


31416 


18850 


6283 


500 


300 


100 1  3000 


Coefficients 


J.A.  ( 

mils) 

6-D 

Theory 

Unstable 

-1.755 
-2. 154 


-1.957 
- 1. 978i 


-1.569 

-2.366 


-1.304 

-2.632 


-0.024 

-3.959 


1.873 

6.399 


3.755- 
10. 393i 

3.751- 
10. 486i 

7. 435- 
34.530i 

7.247- 
34. 867i 

- 1 1 . 870 
-44. 054i 

-12. 144 
-44 . 243i 

-35. 054 
-35. 053i 

-35. 164 
-35. 1 64 i 

TABLE  18.  THEORY  VALIDATION,  ASYMMETRIES, 
CASES  179-189 


TWn  >  WCT  why  ■ 


Figures  22,  23  and  24  illustrate  Cases  157-189.  The  curves^are 
similar  to  those  in  Figures  7,  8  and  9  but  are  displaced  by  the  a 

contribution.  Cases  168,  179,  180  and  181  indicate  that  maximum  effects 
of  the  various  parameters  have  been  achieved  in  other  stable  cases.  The 
effects  of  roll  rate  and  velocity  follow  the  same  trends  as  those  in  Cases 
58-90.  Figure  25  shows  the  effects  of  velocity  for  Cases  157-189.  Cases 
with  U  =  3000  ft/sec  exhibit  the  smallest  dispersicn.  A  sample  trajectory, 
Case  189,  is  shown  in  Figure  26. 

COMPARISON:  11101,  LOW,  VERY  SLOW  ROLL  RATE  THEORIES 

In  Cases  58-189  the  High.  Low,  and  Very  Slow  Roll  Rate  Theories  are 
validated  for  various  initial  conditions  and  parameters.  The  theories 
have  been  applied  for  certain  ranges  in  roll  rate  and  roll  rate  times 
time  (pt) .  The  range  of  pt,  (pt  <_  1.0)  arc  governed  by  the  inherent 
requirements  of  power  series  expansion.  However,  the  ranges  of  p  are 
arbitrary  (to  a  certain  extent)  and  are  based  on  accuracy  of  the  theories 
themselves.  Each  theory  approximates  the  solution  very  well  for  a  certain 
range  of  p  and  then  begins  to  diverge  and  become  inaccurate.  The  range  of 
p  for  which  the  very  slow  roll  rate  theory  is  accurate  is  fairly  well  cut 
and  dried;  p  <  0,  pt  <_  1,0.  For  any  pt  <  1.0  we  must  now  use  the  low  roll 
rate  theory.  “The  question  now  arises,  how  high  a  roll  rate  can  this  theory 
accommodate?  At  what  value  of  p  must  wc  change  to  the  high  roll  rate  theory? 
These  questions  arc  answered  by  a  plot  of  sample  6-11  computations ,  Figure 
27,  and  all  three  theories  extended  beyond  the  limits  used  in  the  previous 
validation.  The  high  roll  rate  theory  is  a  straight  line  going  off  to 
infinity  as  p  goes  to  zero.  Although  the  length  of  the  curve  in  which  it 
is  an  effective  theory  is  short  graphically,  the  range  of  roll  rates  it 
encompasses  is  tremendous.  Figure  28  illustrates  the  effective  limits  of 
each  theory;  that  is,  on  the  spectrum  of  possible  roll  rates  it  shows 

where  each  theory  is  the  most  effective.  The  low  roll  rate  theory  handles 

the  largest  graphical  area  but  only  roll  rates  less  than  100  rad/sec  and 
greater  than  5  rad/sec.  The  upper  limit  of  100  rad/sec  was  chosen  since 

here  the  low  roll  theory  attaches  itself  to  the  0-D  results  while  the  high 

roll  theory  diverges.  The  Lower  limit  of  5  rad/sec  corresponds  to  pt  <_  1.0, 
Figures  27  and  28  depict  Cases  58-68  where  u  =  5000  ft/sec  or  t  -  0.2  sec. 

Therefore  p  =  5  rad/ sec  corresponds  to  pt  =  1.0.  The  very  low  roll  rate 
theory  has  the  smallest  range  but  is  essential  in  predicting  dispersion  as 
the  roll  rate  goes  to  zero.  As  pt  >  1,  the  theory  diverges  as  would  be 
expected  from  a  power  series;  Equation  29.  The  sharp  turn  occurs  at  p  ~  20 
rad/sec  or  pt  ~  4  for  Cases  58-68.  Although  Cases  58-68  were. illustrated 
here,  this  analysis  of  the  effective  limits  of  the  roll  theories  was  found 
to  be  similar  for  all  other  cases.  For  the  u  =  3000  ft/sec  cases  the  low 
roll  theory  limits  were  3.0<p<50  for  u  =  10oB  ft/sec  cases:  1.0<p<25.0. 

PHASE  IV 

To  validate  the  effects  of  gravity  on  dispersion,  a  final  set  of  cases 
were  run  using  the  high  roll  rate  theory,  Equation  24.  Ordinarily,  one 
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Figure  23.  Dispersion:  Phase  III 
Cases  168-178 
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would  think  that  gravity  would  only  introduce  a  constant  terra; 
one  that  could  be  factored  out.  However  integration  of  the 
equations  of  motion  produce  a  gravity  terra  dependent  upon  roll 
rate.  Determination  of  its  validity  and  consequence  is  what 


is  important  here.  So,aQ,  and  aQ  were  set  to  zero  in  order  to 

allow  determination  of  the  effects  due  to  roll  rate  and  veloc¬ 
ity.  The  reduced  governing  equation  becomes: 


and  ot  were  set  to  zero  in  order  to 

o 


_*P*x 

mud 


1000 


No  aerodynamic  asymmetries  were  present  and  the  effects  of 

>,  -j  i 

gravity  were  assumed  independent  of  effects  due  to  So,ao,u  ; 
a  logical  assumption.  Table  19  lists  the  results. 

Table  19  indicates  that  the  effects  due  to  gravity  occur 
largely  in  the  vertical  plane,  as  would  be  expected.  The 
transverse  contribution  is  minimal  but  is  affected  by  both 
velocity  and  roll  rate.  The  vertical  contribution  is  only 
affected  by  velocity.  The  unstable  cases  indicate  maximum 
use  of  magnus  and  thus  maximum  traverse  effects  on  dispersion. 
It  can  be  concluded  from  this  brief  but  through  treatment  that 
gravity  effects  dispersion  only  in  the  vertical  plane  (for  all 
practical  purposes)  and  that  its  contribution  is  constant  with 

velocity.  The  roll  dependent  term,  A,  has  been  shown  to 

mucf 

exist  but  become  negligible  for  the  flechctte.  This  term 
would  possibly  become  important  for  projectile  dispersion  and 
other  missile  applications.  I’rojoetile  motion  with  gravity  is 
typified  by  a  cocking  right  oi  '’o  projectile  in  flight  with  a 
positive  C  but  negative  ;  i,  parameter  A  would  become 

a  ~  n 

negative  and  the  entire  roll  dependent,  term,  positive;  that  is 
cocked  to  the  right,  dispersion  to  the  right.  for  a  finned 
missile  the  opposite  would  occur  due  to  the  agreement  in  sign 
between  . 

■«  J 


TABLE  19.  THEORY  VALIDATION,  GRAVITY 
CASES  190-2Q1 


Initial  Conditions 


PG  uo 
31416 
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31416  ‘ 
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Coefficients 
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SECTION  IV 


FREE  FLIGHT  DATA  ANALYSIS 


In  order  to  analyze  actual  test  firings  as  to  jump  and 
dispersion  and  correlate  them  with  the  validated  theory,  the 
initial  conditions  of  each  test  firing  must  be  obtained  and 
put  into  the  proper  form.  To  obtain  raw  experimental  data, 
test  firings  were  conducted  by  the  U.S.  Army,  Frankford  Arsenal. 
The  configuration  tested  was  the  Producibili ty  Ground  Point 
Flechette,  Figure  29.  The  raw  data  required  was  both  trans¬ 
lational  and  angular;  that  is,  data  was  needed  to  determine 
position  as  a  function  of  time  and  angle  of  attack  of  the 
flechette  as  a  function  of  time.  To  accomplish  this,  Frankford 
Arsenal  devised  the  test  apparatus  shown  in  Figure  30.  The  gun 
barrel  was  mounted  on  a  steel  girder  to  elimiate  "barrel  whip" 
due  to  recoil  and  boresighted  on  a  target  50  meters  down  range. 
At  positions,  1,  3,  5,  7,  9,  and  11  feet  downrange,  dual  high 
speed  cameras  were  placed  to  photograph  the  flechette  as  it 
passed  its  station.  One  camera  was  placed  to  allow  a  top  view 
at  each  station  and  provide  a  means  of  obtaining  swerve  and 
yaw  data.  The  other  camera  allowed  a  side  view  at  each  station 
to  obtain  heave  and  pitch  data.  At  each  station  reference 
marks  oriented  the  flechette  as  to  its  exact  position  downrange. 
This  was  to  allow  for  any  camera  timing  error  and/or  variation 
in  muzzle  velocity. 

From  the  battery  of  testing  firings,  eight  separate  rounds 
were  chosen  at  random  to  be  analyzed.  The  eight  rounds  along 
with  velocity,  roll  rates  and  target  positions  are  given  in 
Table  20. 

Raw  translational  and  angular  data  are  shown  in  Figures  31 
through  46.  The  figures  illustrate  the  position  and  complex 
angle  of  attack  of  the  flechette  for  each  station. 
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gure  29.  Ground  Point  Frechette,  With  and  Without  Sabot 


TABLE  20.  FRANKFORD  TEST  FIRING  DATA 


“o 

Po 

4747 

11,454 

4662 

WEBEGSM 

4642 

14,219 

4662 

13,000 

4758 

13,289 

4753- 

17,354 

Target  at  50  ft 

Y  (ft) 

iZ  (ft) 

0.117 

-0.038 

0.053 

-0.010 

0. 141 

-0.004 

wmsm 

0.099 

0.053 

0.016 

0.084 

wammm 

0.070 

-0.019 

0. 089 

0.059 

16.613 


11,913 


Figure  31 


Raw  Translational  Data 
Ground  Point  Round  4 
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Raw  Translational  Data 
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Raw  Angular  Data 
Ground  Point  Round  16 
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Figure  46.  Raw  Angular  Data 

Ground  Point  Round  19 
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Once  the  raw  data  was  obtained,  it  had  to  be  converted 


into  a  form  such  that  initial  conditions 
extracted  from  it.  To  eventually  arrive 


S„,S  .a  and  c*  could 
o  ’  o  ’  o  _c^ 

at  values  for  and 
o 


be 


S  ,  the  translational  parameters,  the  raw  position  or  transla¬ 
tional  data  had  to  be  approximated  by  equations.  The  raw  data 
was  fitted  to  a  polynomial  equation  of  third  degree  by  a  least 
squares  method.  The  data  in  the  y-direction  was  fit  separately 
from  that  in  the  z-direction  to  distinguish  between  the  swerve 
and  heave  contributions.  With  the  ecuations  obtained,  a  simple 
differentiation  yielded  equations  for  the  _&plocit.iaps  in  the  y 
and  z  directions.  The  initial  conditions  S  and  S  are  now 
readily  obtainable: 


So  <ft)  =  *0  +  iZ° 

SQ  (ft/sec)  =  y 0  +  izQ 

Obtaining  a  and  a  from  the  raw  angular  data  was  more 
”00 

difficult.  The  traditional  way  of  analyzing  any  missile  motion 
with  pitch,  yaw,  and  roll  is  by  a  three-degree-of- freedom  least 
squares  fit  to  the  tricyclic  motion,  Equation  6.  However,  the 
availability  of  only  6  data  points  made  this  technique  impossi¬ 
ble,  so  another,  approximate  method,  had  to  be  employed.  The 
solution  was  to  approximate  the  pitching  and  yawing  motion  to 
one-degree-of- freedom  while  holding  the  roll  rate  constant.  In 
order  to  do  this,  the  3-a  axis  system  had  to  be  rotated  to 
coincide  with  the  more  dominant  angular  mode.  Figure  47  illu¬ 
strates  a  typical  raw  angular  data  plot.  Since  the  angular 
motion  of  the  flechette  tends  to  approximate  an  ellipse,  the 
B-ot  axes  are  rotated  some  angle  y  to  coincide  with  the  major 
and  minor  axes  of  the  ellipse,  as  shown.  The  angular  data  is 
retabulated  for  this  new  axes  system,  3' -a'.  To  fit  the  data 
to  the  one  -  degree  -  of - freedom  equation: 

a  =Kje*'t  cos  (cot  +  6  ) 

only  the  dominant  mode  can  be  considered.  For  example,  in 
Figure  47  the  dominant  mode  occurs  along  the  a'  axis;  therefore, 
only  a'  coordinates  are  utilized  in  the  least  squares  fit, 
corresponding  3'  coordinates  are  ignored.  Table  21  lists  the 
parameters  obtained  for  the  eight  flechette  rounds.  Once  an 
equation  for  a'  is  obtained,  it  represents  one  dimensional 
oscillatory  motion  along  the  a*  axis.  A  simple  differentiating 
of  the  a'  equation  yields  an  equation  for  a'.  The  initial 
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conditions  aQ  and  a Q,  however,  are  complex  whereas  a'  and  a' 

are  only  one  dimensional.  Therefore,  the  rotation  angle  y  is 
taken  into  account  and  the  a'  equation  is  projected  back  into 
the  3,  a  axes  system: 

a=a ' cosy 

a=a ' cosy 

B=a ' siny 

3=a ' siny 


Thus  the  complex  initial  conditions  are  approximated, 
Wloo 


„  .  • 
a  aB  +ia 
o  Mo  o 


Figures  48-63  illustrate  the  fitted  data  both  translational  and 
angular  for  the  eight  rounds.  The  transitional  data  includes 
the  pertinent  equations. 
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igure  48,  Fitted  Translational  Data 
Ground  Point  Round  4 
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Figure  49.  Fitted  Angular  Data  Ground 
Point  Round  4 
95 


T  (FT/SEC)  T  (FEET) 

2.0  1.0  o.o  -l.o  -2.0  0,09  0.08  0.07  0.06  0.0S 


DISTANCE  D0WNRANGE  (FEET) 


oJ 


LlI 

li-nJ 
w0j 

c 


<n 

o 


C3O000fflOO00CTOonaaDOOaDE,t3mQDaQaaa 


Y  -  0.0764  +  0.0519t  -  149.7t2  +  58  94  tv3 
l  m  0.0016  -  1.807t  +  1962.7t?-176737t3 


0  0.0006  0.0012  0.0018  0.0024 

TIME  (SECONDS) 


0.0006 

-  » 


0.0012 


0.0018 

■  ■  I 


0.0024 

—  —  — L..  ....... 


0.0030 


0 . 0030 

_ i- 


Y  •  0.0519  -  299. 4t  +  17682t2 
Z  .  -1.807  +  3925. 4t  -  530211t2 


CJ 

LU 

to_ 

V 


o 

’fMn 


o 

o 


„  n  a  a  Q  a 
aaD 


qbq 


□  a  c 


□  a  Q  a  a° 


a  a  Q 


□  a 


d  oa 


0  Y 

A  z 

0  3  6  9  12 

DISTANCE  D0WNRANGE  (FEET) 
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Figure  51.  Fitted  Angular  Data  Ground 
Point  Round  6 
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Figure  53.  Fitted  Angular  Data 
Ground  Point  Round  7 
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Figure  54.  Fitted  Translational  Data 
Ground  Point  Round  8 
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Figure  55.  Fitted  Angular  Data 
Ground  Point  Round  8 
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Figure  58.  Fitted  Translational  Data 
Ground  Point  Round  16 


104 


.. mi's  S 


(RflO/SEC)  o<  (OEGREES) 

O  -SO  -100  -150  5  4  3  2  10-1-2 


DISTANCE  O0WNRRNGE  (FEET) 

3  6  9 


in-  tr>  ■ 

i  i 

a*.  oi 

'  1  A  A 

cn  cn  r 

i  •  i 


A  *  A 


0  °  a  0  o 


Q 


0.0006  0.0012  0.0010  0.0024  0.0030 

TIME  (SECONDS) 

0.0006  0.0012  0.0016  0.0024  0.0030 


o  I 

ini  a 


0  ®  o 


A  *  A 


Ssl  a  *  *  D 

try  A  A  *  0 


CO  *  A  p 


O-  O 

in  in 

•ft 

o  o 

0‘  o 


o  I  o 
ini  in 


3  6  9 

DISTANCE  O0WNRANGE  (FEET) 

Figure  59.  Fitted  Angular  Data 

Ground  Point  Round  16 


■  WWW* 


Y  (FT/SEC)  Y  (FEET) 

20.0  10.0  0.0  -10.0  -20.0  0.09  0.00  0.07  0.06  0.05 


DISTANCE  O0WNRANGE  (FEET) 

3  6  9  12 

□  Y 

A  z 


15 


o . 

9 

© 


UJ 

LU 

U_fM 

'-'to. 


O 


o  □ 


(9 


Q 


Q 


Q 


*  *  *  A  0 
A  A  .  0 

o$®° 

Q q □ o a  a o o a o a □ 0 0 a Q a 


n 

oJ 


y  «  0.0778  -  26.47t  +  20429t2-  5880767t3+  392393004t'4 
Z  -  0.0006  +  9. 344t  -•  5134t?+  1699768t3-  123239093 t4 


0  0.0006  0.0012  0.0018  0.0024  0.0030 

TIME  (SECONDS) 

0  0.0006  0.0012  0.0018  0.0024  0.0030 


o 

CM 

I 


o 

‘S 

oV 


UJ 

CO 


o 

o 


o 

N° 


o 

CM 


Y  -  -26.47  +  40858t  -  17642301t"+  1569572016t3 
a  Z  -  9.344  -  10268t  +  5099304t2-  492956 364t3  O 


o 


o 


a 


q 


q 


□ 


a 


a 


□ 


o  gj  □ 

- °  °  BP  a  0  0  b  p  a  °  ° 


0 


0 


o 


0 


□  Y 
▲  2 


T"  I  I 


3  6  9  12 

DISTANCE  D0WNRANGE  (FEET) 

Figure  60.  Fitted  Trims  lat  ional  Data 
Ground  Point  Round  17 


15 


106 


in  4 

i 

>  1 

3*  . 

1 

Hi 

CO 

1  9 

1  * 

to? 
UJ  *  ' 

<r>™ 
UJ  1  ' 

LU 

CC7- 

o 

UJ 

sT- 

uj0. 

o° 

S- 

w 

w 

*-«  - 

r*  - 

* 

CM* 

cm" 

CO 

CO' 

*1 

3 '• 

in . 

DISTANCE  OOWNRANGE  (FEET) 

3  6  9  12 

■  I  Im.  ■  A  ■  1,  —  I  I  *  — *  ■ 


15 


0 


0a°0 


0 


0 

0 

□  4 

*  ? 


0 


0 


0 


A  A  A 


0 


0 


0 


0 


0 


0 


A  A 

A  A 

*  A  A  * 


0 


0 


0 


0 


0  0 

0  0 
0  0  0 


0  0.0006  0.0012  0.0018  0.0024  0.0030 

TIME  (SECONDS) 

0  0.0006  0.0012  0.0018  0.0024  0.0030 


0 1  ST  ANCE  DOWNRANGE  (FEET) 

Figure  61.  Fitted  Angular  Data 

Ground  Point  Round  17 


107 


Y  (FT/SEC)  Y  (FEET) 

0.0  10.0  0.0  -10.0  -20.0  0.09  0.08  0.07  0.06  0.05 


[DISTANCE  OQWNRANGE  (FEET) 


6 


12 


15 


0  Y 
A  Z 


o. 

o 


LxJ 
UJ 
Li-  CM 
'-'O. 

o 

rvj 


0000000==>=>oODOEODOOaoSss6 


6. 

a 


cn 

© 


I  Y  .  0.0752  +  7.867t  -  11404t2+  4164544t3--  2952239y6t4 
Z  -  -0.0005  4-  16. 21t  -  6801t?+  1970043t3-  133503655t4 


0  0.0006  0.0012  0.0018  0.002(1  0.0030 

TIME  (SECONDS) 

0  0.0006  0.0012  0.0018  0.0021  0.0030 

■■  -■  .1  I  ■  i  ■■  i  ■!  i  —  .  i  — 1  ■  ■  ■ 


CM 


O 

<N 


1  «  7.867  -  22808t  +  12493632t2-  11808' 5984t3 

2  .  16.21  -  13602t  +  5910129t‘-  534014620t3 


oS 

LU  1 

cn_ 


a 


0 


TS 


-» 


fl  o  II  °  °  °  °  °  °  °  °  a  n 


a- 


nr 


A  A 


A  A 


“  O  _ 

0  D 

***“****»di*4l 

□  *  A  A 


0 


0 


0 


O 

CN 


□  Y 
A  Z 


a 


3  6  9  12 

DISTANCE  DOWNRANGE  (FEET) 


□ 


'  I  ■  "  ^ 

15 


liiMiif  hi.  Pitted  Translational  Data 
Ground  Point  Round  ID 


l  ns 


DISTRNCE  OOWNRRNGE  (FEET) 

3  6  9 


0  ct 

A  f 


Am®G,0°D  4 
A  O  m 

O 

AO  Q 

o 


TABLE  21.  AERODYNAMIC  PARAMETERS  FROM 
LEAST  SQUARES  FIT 


R 

0 

U 

N 

D 

Ki 

(degrees) 

A 

(rad/sec) 

& 

(rad/sec) 

6 

(rad) 

4 

5.01 

-49.48 

1921.3 

-1.29 

6 

3.64 

68.24 

2079.8 

-1.21 

7 

-2.78 

46.48 

1871.4 

-1.26 

8 

-4.02 

-203.  19 

2267.6 

-1.21 

14 

6.09 

-126.37 

2042.0 

-1.23 

16 

5.84 

-174.7 

'5211.7 

-1.  18 

17 

-4.81 

8.53 

23 14. 5 

-1.02 
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1889.9 
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a=K^e  cos(wt  +  6) 
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SECT  I  ON  V 

DISPERSION  ANALYSIS 


FREE  FLIGHT  VERSUS  THEORY 


Once  the  initial  conditions  arc  determined  as  in  the  pre¬ 
vious  section,  they  are  applied  to  the  theory  and  compared  to 
the  dispersion  of  each  test  fired  round.  To  utilize  the  theory, 
the  fitted  data  must  be  chosen  for  a  given  time;  that  is, 


So,®o,clo»®o  must  selected  for  one  given  point  in  time-posi¬ 
tion  downrange.  Since  the  question  of  what  point  in  time 
the  initial  conditions  occur,  3  sets  of  initial  conditions  were 


chosen  to  correspond  with  positions  1,  3,  5  feet  downrange. 

This  span  of  position  downrange  may  or  may  not  be  sufficient  to 
include  the  actual  time  corresponding  to  the  initial  conditions 
for  each  round.  The  following  analysis  will  determine  each 
round's  effective  time  for  its  initial  conditions. 

For  each  set  of  initial  conditions,  theory  and  6-D  compu¬ 
tations  were  done  and  compared  to  target  data  for  the  Frankford 
test  firings.  The  results  are  tabulated  in  Table  22  in  mils 
and  plotted  in  Figures  64-71  in  feet;  deviation  from  the  time 
of  fire  at  50  feet  downrange.  The  relationship  between  the 
deviations  in  feet  and  mils  at  50  feet  downrange  is: 

17a.  (mils)  *  JLlftl  (1000) 
x 


or  J.  A.  (mils)  =  (20)  S  (ft) 

To  accurately  and  concisely  analyze  the  complex  and  large 
amount  of  data  in  Table  22,  the  positions  downrange  in  which 
the  initial  conditions  were  selected  must  be  simultaneously 
analyzed  with  the  dispersion  results  at  50  feet  downrange.  The 
problem  in  choosing  initial  conditions  is  where  they  should  be 
taken;  at  what  point  downrange.  Normally,  one  would  think  that 
the  initial  conditions  would  occur  immediately  after  leaving 
the  gun  barrel.  However,  the  flcchette  being  a  finned  body 
needs  a  sabot  configuration  to  guide  it  down  the  barrel,  Figure 
29.  The  sabot  causes  the  initial  condition  location  problem 
since  the  sabot  must  separate  from  the  flechette  outside  of  the 
gun  barrel.  The  exact  time  and  place  where  this  occurs  is  not 
constant;  varying  from  round  to  round.  Not  only  does  the  sabot 
separate  from  the  flcchette  instantaneously  different  every 
time,  the  sabot  may  not  separate  cleanly  or  the  same  way  every 
time.  Interference  with  the  fins  after  sabot  separation  can 
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TABLE  22.  DISPERSION  ANALYSIS  RESULTS 


INITIAL  CONDITIONS  TAKEN  AT 
1,3,5  FEET  OOWNRANGE 


■  TEST  FIRING  CD  THEORY 


Dispersion:  Ground  Point-Round  4  Tost  Firing 

Versus  Theory  at  50  Feet  Downrange 
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-FEET 


INITIAL  CONDITIONS  TAKEN  AT 
1,3,5  FEET  OOMNhAHGE 

■  TEST  FIRING  0  THEORY 


INITIAL  CONDITIONS  TAKEN  AT 
1,3,5  FEET  OOWNRANGE 


■  TEST  FIRING  □  THEORY 


Figure  68.  Dispersion:  Ground  Point  Round  14  Test 

Firing  Versus  Theory,  at  50  Feet  Downrange 
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INITIAL  CONDITIONS  TAKEN  AT 
1.3.5  FEET  OOWNRANGE 

■  TEST  FIRING  □  THEORY 


Dispersion:  Ground  Point  Round  17  Test 

Firing  Versus  Theory  at  50  loot  Downrangc 


Z-FEET 


Y-FEET 


-0.2  0.1  0  0.2  0.M  0.6 


Figure  71.  Dispersion:  Ground  Point  Round  19  Test 

Firing  Versus  Theory  at  50  Feet  Downrange 
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cause  d i s tu r banco s  to  the  flechette  and  alter  the  initial  con¬ 
ditions.  In  addition,  asymmetric  sabot  separation  can  influ¬ 
ence  the  initial  conditions.  Therefore,  the  positions  down- 
range  of  1,  3,  5,  and  7  feet  were  photographed  to  show  the 
flight  transition  sequence.  Figures  72-79  illustrate  the 
flight  transition  sequence  for  the  8  flechette  test  rounds. 

In  every  sequence  the  sabot  begins  to  separate,  in  varying 
degrees,  1  foot  downrange.  At  3  feet  downrange,  the  sabot  is 
nearly  completely  separated,  but  in  some  cases  the  sabot  par¬ 
ticles  pose  interference  problems  with  the  fins.  By  5  and  7 
feet  downrange  the  sabot  has  completely  separated  and  the  fle¬ 
chette  is  in  free  flight.  The  correspondence  between  the 
flight  transit  ion  sequence  and  dispersion  results  can  be  seen 
in  each  individual  round.  Figure  64  indicates  that  the  initial 
conditions  for  round  4  occur  somewhere  between  1  and  3  feet 
downrange  judging  by  the  dispersion  of  the  actual  tested  round. 
Figure  72  verifies  this  fact  in  that  the  sabot  has  separated 
from  the  flechette  between  1  and  3  feet  downrange.  The  y-coordi- 
nate  in  the  dispersion  vector  does  not  accurately  agree  with 
the  theory  for  this  case.  However,  besides  computational  error 
other  physical  factors  can  influence  dispersion.  Contributions 
by  fin  asymmetries  and  other  configurational  as/mmetries  can  be 
important  but  are  unable  to  be  detected  or  accounted  for. 
Throughout  this  analysis  this  must  be  kept  in  mind  to  partially 
account  for  any  discrepancy  between  the  actual  test  firing  and 
the  theory  and  0-D  computations.  Figure  65  indicates  the  ini¬ 
tial  conditions  for  round  6  occur  between  3  and  5  feet  down- 
range.  Figure  73  verifies  this  choice  showing  separation 
occurring  around  3  feet  but  with  sabot  particles  very  close  to 
the  fins  causing  possible  interference  and  delaying  the  initial 
conditions  location.  The  initial  conditions  location  for  round 
7  is  difficult  to  accurately  choose  since  the  y-coordinate  does 
not  accurately  agree,  Figure  66.  It  is  safe  to  say  that  the 
initial  conditions  occur  sometime  around  3  feet  and  Figure  74 
verifies  this  choice.  The  z-coordinate  for  round  8  is  not  as 
accurate  as  would  be  desired,  Figure  7,  but  the  y-coordinate 
indicates  initial  conditions  occurring  between  3  and  5  feet 
dowhrange.  Figure  75  agrees  with  this  choice  indicated  inter¬ 
ference  witli  the  fins  at  3  feet  delaying  the  initial  conditions. 
Initial  conditions  for  round  14  are  chosen  between  3  and  5  feet 
downrange,  Figure  68.  Figure  76  indicates  possible  fin  inter¬ 
ference  tending  to  verify  the  choice.  Figures  69  and  77  indi¬ 
cate  and  verify  the  choice  of  initial  conditions  in  the  immedi¬ 
ate  vicinity  of  3  feet  downrange  for  round  16.  Possible  fin 
interference  at  3  Feet  downrange,  Figure  78,  round  17,  verifies 
a  choice  of  initial  conditions  between  3  and  5  feet,  Figure  70. 

A  similar  situation  occurs  for  round  19  in  Figures  71  and  79. 

It  is  often  difficult  to  choose  initial  condition  positions 
accurately  due  to  slight  discrepancies  between  theory  and  test 
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Figure  76.  Flight  Transition  Sequence  Round  14 


TCP  VIEW 


firings.  However,  the  discrepancies  are  of  the  order  0.05  feet, 
which  shows  up  large  in  Figures  64-71  due  to  the  scale  chosen, 
but  is  within  the  error  expected  from  the  validation  of  theory 
section . 

The  influence  of  sabot  separation  can  be  readily  seen  by 
inspection  of  Figures  72-70,  1  and  3  feet  downrange.  In  every 
case,  the  flechette  and  sabot  are  at  nearly  a  zero  angle  of 
attack  at  1  foot,  but  has  changed  angle  of  attack  noticeably  by 
3  feet  downrange.  This  would  indicate  that  fin  interference  or 
asymmetric  sabot  separation  is  causing  the  noticeable  effect. 

It  can  be  concluded  that  dispersion  is  dependent  upon  the  ini¬ 
tial.  conditions  that  the  initial  conditions  are  a  function  of 
sabot  separation  and  that  the  theo-y  ca:  predict  what  the  ini¬ 
tial  conditions  are  and  where  they  occur. 

DISPERSION  THEORY  VERS.S  FIRST  MAXIMUM  YAW  HYPOTHESIS 

A  popular  theory  to  predict  the  dispersion  of  flechette  is 
the  First  Maximum  Yaw  Hypothesis.  This  theory  relates  the  dis¬ 
persion  magnitude  co  the  first  maximum  yaw  magnitude  by  a  nearly 
linear^relationship .  Other  initial  conditions  such  as  angular 

rate,  aQ,  and  translational  velocity,  SQ  are  said  not  to  effect 

dispersion.  To  disprove  this  theory  and  strengthen  the  posi¬ 
tion  of  the  theory  ascribing  to  dispersion  due  to  initial  condi¬ 
tions  the  First  Maximum  Yaw  theory  was  applied  to 

Frankford  Arsenal  data.  Figure  80  shows  a  plot  of  disperson 
magnitude  versus  first  maximum  yaw  magnitude.  Clearly  no  linear 
relationship  exists  between  dispersion  and  first  maximum  yaw. 

In  fact,  the  plotted  data  resembles  a  random  shotgun  blast. 
Figures  81,  82,  and  83  employ  the  theory  to  the  first  maximum 
yaw  hypothesis.  Again  the  plot  substantiates  the  findings  of 
Figure  80.  The  uisproval  of  the  first  maximum  yaw  hypothesis 
comes  as  no  surprise  since  the  dispersion  theory  contradicts  it 
and  the  6-D  computations,  which  integrate  the  actual  equations 
of  motion,  validated  the  dispersion  theory.  Therefore,  disper¬ 
sion  could  never  accurately  be  predicted  by  a  theory  involving 
only  first  maximum  yaw. 

_x  _v 

The  influence  of  initial  conditions,  S  ,a  ,  and  a  and 

’  o  o '  o 

dispersion  for  the  actual  test  firings  are  expected  to  be 
different  from  that  in  the  validation  of  theory  section  because 

of  the  different  ranges  in  the  initial  conditions.  For  example, 
• 

S  only  in  the  validation  sect  ion  was  (100  +  lOOi)  ft/sec.  In 

O  -—Vi 

the  actual  test  firings,  only  ranged  as  high  as  20.0  ft/sec. 
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Of  course,  the  large  value  was  only  to  validate  the  theory. 

■  »■  ^ 

However,  SQ  is  a  very  important  contributor  to  dispersion.  In 

the  reduced  equation  24,  employed  to  calculate  the  theory  column 
in  Table  22, 


J.  A.(mils)=  1000 

for  round  4,  1  foot  downrange, 


O  S0  Iy  (-+• 

—  +  ~  A  a  -  a  ■ 

u  mud  V  o  o 


iplx, 
I, 


1000 


u 


=  Cl  .  10  0  7  +  1  .91  34  i.)  mils 


where  as, 


J.  A.  *  (1.329  •  1.302i)  mils 


Since  this  is  typical  of  the  8  rounds  tested,  S  has  a  profound 
effect  on  dispersion  for  these  rounds. 

Similarly,  for  this  particular  case, 


1000 


Jo 


(O.OlbO  +  0.0004  i)  mils 


ipIvA 

1000  a  JUL-  =  (-0.01437  +  0.00214i)  mils 
o  mud 


.  j 

"100°  mud  =  ('°-206075  •  1.3836721)  mils 

Obviously,  s  and  aQ  are  by  far  the  greatest  contributors  to 

dispersion  for  this  case.  Inspection  of  all  the  other  31  cases 
in  Table  22  agrees  with  this  general  pattern.  "5  can  be  nearly 

eliminated,  of  course,  by  accurate  setup  of  the  test  equipment 
so  that  the  gun  barrel  is  set  exactly  at  coordinates  (0,0). 

Any  IT  then  would  occur  from  displacement  due  to  the  blast. 

figure  84  illustrates  the  dependence  of  the  Jump  Angle,  and  hence 
dispersion,  upon  angular  rate  and  angle  of  attack. 


Although  SQ  and  aQ  contribute  the  most  to  the  dump  Angle, 

the  combination  of  SQ  and  «0  can  have  a  noticeable  influence. 

From  the  test  firings,  S  was  found  to  have  a  negligible  effect 

on  dispersion.  Therefore,  it  is  neglected  in  Figure  84  to 
simplify  the  It  is  evident  from  Figure  84  that  -various 

combinations  of  aQ  aQ  yield  zero  dispersion.  It  is  possible 
that  iargc  values  of  a  and  aQ  can  combine  to  yield  zero  disper¬ 
sion;  an  impossibility  with  the  first  maximum  yaw  hypothesis. 

If  a  and  aQ  are  able  to  balance  to  give  zero  dispersion,  then 
this  idea  cun  be  expanded  to  include  the  entire  equation. 

The  governing  equation  used  throughout  this  dispersion 
analysis  section  is: 

— *-  S0  Iy  *Plx\  ig  /  x  \ 

J.A.-  1000  —  +T  “  mud  A(^o~ao  ^  )+  2  \^2~) 


Eliminating  the  constant  gravity  term, 


~~  ■  hm  —y  ^ 

JTa.  =  1000  ^  +  ?2.  -  Ll. -a  fill  -  cu  ^ 

J  x  u  mud  l  o  uo  i 


Setting  rt.  to  zero,  the  idea  behind  Figure  84  is  expanded  to 


include  S  .  S  . 

o  o 


fo  +  h  =  a  (~6t 

x  u  mud  ao  0  Lr  j 


rearranging 


Xf) +  X  = 


<“o  Vao  ‘Pit' 


A  dimensional  analysis  of  the  equations  finds  that  both  sides 
have  units  of  momentum.  Going  one  step  farther  it  can  be  said 
that  to  obtain  zero  dispersion:  initial  transverse  momentum9 
initial  momentums  that  causes  dispersion.  The  size  of  initial 
conditions  can  be  huge,  Figure  84,  but  if  they  can  combine  to 
balance,  zero  dispersion  results.  The  way  the  initial  condi¬ 
tions  combine,  determine  the  magnitude  of  the  imbalance  or 


dispersion.  It  should  be  noted  that  this  dispersion  discussed 
is  round  to  round  dispersion  and  that  the  inconsistency  of  the 
momentum  imbalance  from  round  to  round  causes  a  dispersion 
pattern  (a  set  of  rounds) .  The  next  section  will  highlight 
this  principle  in  the  evaluation  of  physical  factors  affecting 
dispers ion . 
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SECTION  VI 


PHYSICAL  EVALUATION  OP  DISPERSION 


Initial  momentum  imbalance  has  been  shown  to  cause  disper¬ 
sion.  Initial  conditions  determine  the  magnitude  of  the  imbal¬ 
ance.  What  causes  these  initial  conditions  to  occur  is  the 
subject  of  this  final  section.  Initial  conditions  occur  some¬ 
where  between  0  and  5  feet  downrange  to  different  degrees  of 
magnitude  due  to  various  conditions.  These  conditions  are: 

1.  Fin  or  body  asymmetry 

2.  In-bore  mal - al ignment 

3.  Asymmetric  blast 

4.  Asymmetric  sabot  separation 

5.  Sabot-fin  interference 


6.  Pin  or  body  damage 

Fin  or  body  asymmetries  can  cause  dispersion  magnitudes  to  range 
as  much  or  greater  than  those  in  the  Validation  of  Theory  section 
for  aerodynamic  asymmetries.  These  asymmetries  can  be  overcanted 
or  bent  fins,  damaged  nose  cone,  or  even  body  deformities.  Fig¬ 
ure  85  which  shows  in-bore  mal - al igmnent  also  shows  a  slightly 
bent  body,  concave  downward.  In-bore  mal -alignment  can  be 
attributed  to  warping  and/or  the  entire  flechette  at  some  angle 
of  attack.  Clearly,  if  this  flechette  were  fired,  the  in-boro 
angle  of  attack  would  produce  an  ‘o’  outside  the  gun  barrel  even 

before  sabot  separation.  With  the  flechette  at  some  angle  of 

*-X  __k. 

attack,  the  blast  can  cause  a  large  a  and  an  and  S  .  The 

blast  itself  is  a  chief  catalyst  in  causing  the  initial  condi¬ 
tions.  An  asymmetric  blast  can  indeed  impart  influence  on  the 
initial  conditions,  but  a  symmetric  blast  can  also.  Given  an 
initial  angle  of  attack  due  to  some  disturbances  the  symmetric 

blast  can  cause  significant  o'  ,ct„,S  and  Figure  86  shows 

a  typical  blast  region  with  the  flechette  outlined  in  the  pic¬ 
ture.  The  momentum  principle  discussed  in  the  previous  section 
goes  hand-in-hand  with  this  blast  region.  It  is  here  that  the 
transverse  and  angular-momentum  is  imparted  to  the  flechette. 
Figure  87  illustrates  a  typical  flechette  in  the  blast  region. 
Coming  out  of  the  barrel  at  some  angle  of  attack,  the  blast 
catches  the  flechette  and  induces  some  angular  rate.  At  the 
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same  time,  the  flechette  is  translated  laterally  giving  an  S 
m ,  o 

and  Sp.  If  these  contributions  cancel  each  other  out;  that  is, 
if  initial  transverse  momentum  equals  initial  angular  momentum 
then  the  dispersion  is  zero.  If  they  do  not  cancel,  dispersion 
results.  The  sketch  is  highly  simplified  in  that  the  blast 
itself  is  all-engulfing  as  in  Figure  86.  Of  course,  the  transi 
tion  sequence  of  sabot  separation,  fin  interference,  and  possi¬ 
ble  fin  damage  must  not  be  forgotten.  The  transition  sequence 
occurs  in  the  blast  region,  however,  and  is  not  considered 
separate  from  the  blast.  When  separation  occurs,  the  sabot 
particles  are  apt  to  interfere  with  the  fin  section  and  cause 
possib1^  damage.  Once  the  sabot  has  separated  and  cleared  the 
fins,  the  blast  has  had  its  greatest  effect  and  the  initial  con 
di tiers  can  be  determined.  After  the  flechette  has  moved  down- 
i ange ,  it  assumes  a  supersonic  free  flight,  Figure  88. 


Flechette  In-Bore  Position 


Figure  86.  Typical  Flechette  Blast  Region 


Figure  87.  Muzzle  Blast  Effects 
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SECTION  VII 


CONCLUSIONS 


A  complete  Jump  and  Dispersion  Theory  has  been  developed 
for  free  flight  vehicles.  Three  governing  equations  have  been 
determined  to  accommodate  high,  low,  and  very  low  roll  rates. 

The  theories  were  found  to  be  accurate  with  six-degree-of - free  - 
dom  numerical  computations  of  the  equations  of  motion  and  there¬ 
fore  reliably  predict  the  jump  and  dispersion  of  l'lechettes. 

The  theory  validation  included  201  case  runs  in  four  phases. 

The  first  phase  validated  the  theory  with  respect  to  restoring 
and  damping  moments.  The  effect  of  these  moments  on  dispersion 
was  found  to  depend  on  the  initial  conditions.  The  second  phase 
validated  the  theory  with  respect  to  Magnus  forces  and  moments. 
The  effect  of  Magnus  was  found  to  be  very  small  and  not  of  any 
consequence  unless  the  total  dispersion  of  any  given  round  was 
of  the  same  order  of  magnitude  as  the  Magnus  effect.  Phase 
three  validates  the  theory  with  respect  to  aerodynamic  asymmet¬ 
ries  and  roll  rate.  All  three  theories  were  validated  in  this 
phase  and  found  quite  accurate  considering  the  large  dispersions 
encountered.  Aerodynamic  asymmetries  causing  a  trim  angle  of 
1  degree  had  little  effect  on  the  dispersion  of  flechettes. 
Slower  rolling  bodies  were  shown  to  have,  in  general,  increas¬ 
ingly  larger  dispersion  values  as  roll  rate  decreased.  It  can 
be  concluded  that  for  free  flight  vehicles  prone  to  aerodyna¬ 
mic  asymmetries  and  fin  damage,  a  high  roll  rate  is  essential 
to  lower  dispersion  and  increase  accuracy.  The  fourth  phase 
validates  the  theory  with  respect  to  gravity.  The  theory 
indicates  a  lateral  contribution  to  dispersion  from  gravity  in 
addition  to  the  obvious  vertical  contribution.  For  the  fle- 
chette,  the  lateral  contribution  was  found  to  be  minimal  and 
was  neglected  in  this  analysis. 

Free  flight  data  was  obtained  from  Frankfnrd  Arsenal  to 
correlate  with  the  theory.  Angular  and  translational  data  was 
fitted  and  put  into  initial  condition  form.  The  initial  condi¬ 
tion  data  was  applied  to  the  theory  and  compared  to  target  data 
for  the  rounds  tested.  The  theory  was  found  to  agree  favorably 
in  magnitude  with  the  test  firings.  As  a  result,  the  method 
used  to  analyze  the  data  can  be  considered  a  valid  method. 
Photographs  of  the  test  firings  were  taken  to  include  the  flight 
transition  sequence  in  the  blast  region.  The  pictures  further 
verify  the  analysis  method  of  the  initial  conditions  by  allow¬ 
ing  agreement  between  the  chosen  initial  conditions  and  the 
position  downrange  where  they  were  selected. 

The  evaluation  of  the  free  flight  dispersion  against  the 
theory  also  disproves  the  First  Maximum  Yaw  hypothesis.  A  plot 
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of  jump  angle  versus  first  maximum  yaw  of  actual  test  data 
produced  a  shotgun  blast  pattern  with  no  relationship  evident 
between  dispersion  and  first  maximum  yaw.  In  addition,  a  plot 
of  jump  angle  versus  angular  rate  for  various  initial  angles 
of  attack  indicates  an  infinite  amount  of  combinations  of 
initial  conditions  to  yield  a  given  jump  angle.  Thus,  zero 
dispersion  has  an  infinite  set  of  possible  initial  conditions. 
It  was  found  for  zero  dispersion  that  a  unique  physical  condi¬ 
tion  holds:  to  obtain  zero  dispersion,  initial  transverse 
momentum  =  initial  angular  momentum.  These  momenta  are 
imparted  to  the  flechette  in  the  blast  region  where  the  body 
and  especially  the  fins  are  subject  to  disturbances.  Momentum 
imbalance  is  the  reason  dispersion  occurs.  The  initial  condi¬ 
tions  only  determine  the  magnitude  of  imbalance  or  dispersion. 
This  dispersion  is  round  to  round  dispersion.  Inconsistency 
in  the  imbalance  results  in  a  dispersion  pattern.  The  initial 
conditions  were  found  not  to  occur  until  after  the  sabot  sepa¬ 
ration  and  the  blast  has  had  its  greatest  effect.  The  factors 
causing  the  existence  of  initial  conditions  were  found  to  be 
not  only  the  blast  and  sabot  separation  sequence,  but  also  fin 
and  body  asymmetries  and  bore  mal-alignment .  In  order  to 
decrease  dispersion,  these  physical  factors  causing  initial 
conditions  must  be  kept  at  a  minimum.  The  most  important 
aspect  would  be  to  protect  the  fins  from  asymmetries,  damage, 
and  interference  from  the  separating  sabot.  These  initial 
conditions  can  never  realistically  be  eliminated  but  if  kept 
minimal,  dispersion  is  reduced. 
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APPENDIX 


This  appendix  contains  mass  parameters  and  stability 
coefficients  for  the  Ground  Point  Flechette.  Table  A-l  lists 
values  for  mass,  dismeter,  axial  and  transverse  moments  of 
inertia.  Figures  A-l  through  A-8  present  stability  coeffi¬ 
cients  used  in  this  analysis  versus  Mach  number.  C  ,  C.,  , 

a  a 

+  Cw.  were  provided  by  Frankford  Arsenal.  C  ,  CM  , 

Mq  p3  p3 

^YE’  ^ZE’  ^ME’  ^NE  were  nominal  values  of  the  coefficients 

following  the  same  trends  of  Cz  and  for  Mach  number. 

a  1  a 

and  were  verified  in  the  University  of  Notre  Dame 

a  q  a 

supersonic  wind  tunnel  (Reference  16)  . 

TABLE  A-l 

FLECHETTE  PARAMETERS 


mass  =  0.000046  slugs 

diameter  =  0.006  feet 

I  =  0.000000000217  slugs-ft2 

I  =  0. 000000036421  slugs-ft2 
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NORMAL  FORCE  COEFFICIENT 


2.00  3.00  14.00  5.00 

RCH  NUMBER 
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R3YMMETRY  FORCE  COEFFICIENT 


Figure  A-6.  Cyg  ,  C^p  Versus  Mach  Number 
Producib i 1 ity  Ground  Point 
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